Relativistic effects in galaxy clustering in a parametrized post-Friedmann universe 
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We explore the signatures of quintessence and modified gravity theories in the relativistic descrip- 
tion of galaxy clustering within a parametrized post-Friedmann framework. For this purpose, we de- 
velop a calibration method to consistently account for horizon-scale effects in the linear parametrized 
Post-Friedmann perturbations of minimally- and non-minimally-coupled scalar-tensor theories and 
test it against the full model-specific fluctuations. We further study the relativistic effects in galaxy 
clustering for the normal and self-accelerating branches of the Dvali-Gabadadze-Porrati braneworld 
model as well as for phenomenological modifications of gravity. We quantify the impact of modified 
gravity and dark energy models on galaxy clustering by computing the velocity-to-matter density 
ratio T , the velocity contribution 1Z, and the potential contribution V and give an estimate of their 
detectability in future galaxy surveys. Our results show that, in general, the relativistic correction 
contains additional information on gravity and dark energy, which needs to be taken into account 
in consistent horizon-scale tests of departures from ACDM using the galaxy-density field. 
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I. INTRODUCTION 

Modifications of gravity and dark energy models can 
serve as an alternative explanation to the cosmological 
constant causing the late-time accelerated expansion of 
our universe. There exist a large number of models that 
yield cosmic acceleration and it is therefore of great in- 
terest to cultivate a formalism within which such models 
can systematically be explored while providing a con- 
sistent description of their small- and large-scale struc- 
ture. Whereas the background evolution can be nearly 
identical for all models, their perturbations generally dif- 
fer. In 'analogy' to the parametrized post-Newtonian 
formalism on local scales [l|, in the last few years, a 
parametrized post-Friedmann (PPF) framework [95| has 
been developed for the generalized computation of the 
linear perturbations around the spatially ho mog eneous 
and isotropic background of our universe |2-KL8j. The 
PPF formalism consistently unifies this computation for 
a large class of modified gravity and dark energy models 
on all linear scales, i.e., if correctly accounting for near- 
horizon modifications, and can be used to consistently 
test these models with structures from the largest ob- 
servable scales. Hence, it provides the proper generalized 
framework to study effects on the relativistic description 
of the galaxy-density field for modified gravity and dark 
energy cosmologies. 

Among the most extensively studied models of modi- 
fied gravity are scalar-tensor theories. Such models ap- 
pear naturally in attempts to unify general relativity 
(GR) with the standard model interactions. In these 
models a scalar field appears in addition to the gravi- 
tational tensor field. This fifth element or quintessence 
field [HI may not couple directly to the matter fields 
but can source the cosmic acceleration observed today. 
Thus, such minimally-coupled scalar-tensor models are 
natural candidates for dark energy. The scalar field 



may, however, couple to the matter fields non-minimally, 
yielding modifications of the gravitational interactions. 
Such extended or generalized quintessence models have 
intriguing properties and may, e.g., act as the source 
of a phantom energy, where the equation of state as- 
sociated with the scalar field becomes of the form p < 
—p. Non-minimally-coupled scalar-tensor theories are 
strongly limited by local measurements (see, e.g. [l[). 
However, nonlinear mechanisms may appear that shield 
the extra force in such high-density environments and al- 
leviate the constraints (see, e.g., [Io|). In the light of 
known (see also [2l|) and possibly yet unknown suppres- 
sion effects, it is important to independently test gravity 
on large scales in order to confirm the extrapolation of 
GR from the well-tested local region to the largest scales 
of our universe. 

With the advent of large-scale galaxy redshift surveys 
in the past few decades, galaxy clustering has become 
one of the most powerful measures of the large-scale in- 
homogeneity and probes GR on scales untested before. 
In particular, in [22} , measurements of galaxy clustering 
in rcdshift-spacc from the Baryonic Oscillation Spectro- 
scopic Survey (BOSS) have been used to constrain devi- 
ations from GR parametrized by the logarithmic growth 
rate of structure / = Q^. Furthermore, constraints on 
the redshift-evolution of / are inferred in (23j utilizing 
data from the WiggleZ Dark Energy Survey. All of their 
measurements are well described by a ACDM universe 
with GR. 

It is important to keep in mind, however, that on larger 
scales, close to the horizon, where modified gravity or 
dark energy may drive the cosmic acceleration today, 
the standard Newtonian description of galaxy cluster- 
ing becomes inaccurate as the relativistic effects become 
substantial. Measurements on these large scales may 
be misinterpreted as the breakdown of GR on horizon 
scales. Moreover, certain classes of models such as Dvali- 
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Gabadadze-Porrati (DGP) [24[ braneworld gravity pos- 
sess a characteristic scale at which the graviton starts to 
propagate into the codimensions of the bulk embedding 
our 4D brane universe and features a transition from the 
effective gravity model in 4D to its full nature of modified 
gravity. Therefore, testing GR against those models re- 
quires a proper covariant description of galaxy clustering 
and its measurements beyond the transition scale. 

In the past few years, the relativistic description of 
galaxy clustering has been developed |25|, [26[ to take 
into account the relativistic effects in galaxy clustering at 
large scales. It is shown (2514321 ] that on large scales, the 
signature of galaxy clustering in GR is significantly dif- 
ferent from its Newtonian description and that the rela- 
tivistic effects provide a new horizon-scale test of GR [33| . 
The relativistic formula differs from the Newtonian de- 
scription because our measurements of galaxy clustering 
are constructed by observing photons and they are di- 
rectly affected by the two metric potentials. Compared 
to ACDM, in modified gravity and dark energy models, 
the metric potentials and the velocity respond differently 
to the same matter distribution. This difference charac- 
terizing the modifications to ACDM may be prominent 
on very large scales, providing a unique opportunity to 
probe gravity and dark energy. 

In this paper, we explore the observable signatures of 
modified gravity and dynamical dark energy in galaxy 
clustering. We adopt the rather phenomenological but 
simple and efficient PPF formalism of @ (cf. Qjj : for a 
unified and consistent description of the associated lin- 
ear perturbations. For generalized scalar-tensor theories, 
horizon-scale PPF fits within the formalism of m have 
not been developed previously (see, however, [1J| f° r a 
description of scalar-tensor model perturbations in a dif- 
ferent PPF approach). Therefore, in order to consistently 
compute the near-horizon perturbations, we introduce a 
calibration method for these models. We test our fits 
against the full model-specific linear perturbations and 
compare our results to existing PPF fits in the limit of 
f(R) gravity, corresponding to a specific case of scalar- 
tensor gravity. We then generalize the relativistic formula 
of galaxy clustering for the description of a PPF uni- 
verse and explore the signature of scalar-tensor gravity 
and quintessence models in galaxy clustering, estimating 
their detection significance in future galaxy surveys. In 
addition to the scalar-tensor models, we give predictions 
for the self- accelerating and normal branch DGP models 
as well as for phenomenological modifications of gravity, 
all described through the PPF formalism. 

The organization of the paper is as follows. In mi Al 
we first review the most important aspects of minimally- 
and non-minimally-coupled scalar-tensor gravity theo- 
ries, present specific examples of such models, and study 
their background evolutions. We devote ^11 Bl to the 
DGP braneworld model and ^11 CI to the phenomenolog- 
ical modifications of gravity. In mill we review the PPF 
formalism of @ describing the linear perturbations on the 
cosmological background within these models. In ijlllBl 



wc describe our calibration method for computing the 
near-horizon perturbations in scalar-tensor models and 
discuss its performance in comparison to the full pertur- 
bation analysis. In mil CI and mil Dl we present the PPF 
fitting functions and parameters for f(R) gravity and 
DGP [(| Hi| , respectively. The phenomenological modifi- 
cations of gravity are defined through the PPF formalism 
in mil El In mVl wc discuss the application of the PPF 
formalism to determine the signature of modified grav- 
ity and dark energy models in galaxy clustering. We 
quantify their impact on galaxy clustering in QIV Al and 
present approximate estimates of their detectability in 
future galaxy surveys in i flVBI We conclude in ffV] and 
refer to Appendix |A1 and 151 for a review on the linear per- 
turbation theory of scalar-tensor models and the details 
on the codes employed for the numerical computations. 

Throughout the paper, we assume a spatially homo- 
geneous and isotropic universe. Hence, the background 
metric is defined by the Friedmann-Lemaitre- Robertson- 
Walker (FLRW) line element 

ds 2 = -dt 2 + a 2 (i)dx 2 , (1) 

where we further assume spatial flatness K = 0. 
We adopt h = 0.73 and f2 m = 0.24 based on the 
ACDM cosmology inferred from the Wilkinson Mi- 
crowave Anisotropy Probe [35| in all modified gravity and 
dark energy models throughout the paper to highlight 
deviations from a ACDM universe. 



II. DARK ENERGY AND MODIFIED GRAVITY 

In order to assess the effects of modified gravity and 
dark energy models on the relativistic description of 
galaxy clustering, wc specialize to a number of well- 
studied quintessence and modified gravity models along 
with phenomenological modifications defined through the 
PPF formalism. We summarize the models studied in 
this paper in Table U (see £|IVB[) . In this section, we 
shall give their background equations with a brief model 
description, while their perturbation equations are de- 
scribed in 31111 We devote ^11 Al to the minimally- and 
non-minimally-coupled scalar-tensor theory. In ^11 Bl 
and mi CI we discuss the DGP braneworld scenario and 
phenomenological modifications of gravity, respectively. 

A. Scalar-tensor gravity 

We consider scalar-tensor theories in which the modi- 
fied Einstcin-Hilbcrt action in the Jordan frame is of the 
form 

S = ^ J d i x^[F{ V )R-Z{ V )d^d^-2U( V )\ 

}, (2) 

where S m is the matter action with matter fields i/j m , 
k' 2 = 8irG with the bare gravitational coupling G, and 
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we have set the speed of light c = 1. Here, F(<p) describes 
the coupling of the scalar field tp to the metric, Z(tp) is 
the kinetic coupling, and U(tp) is the scalar-field poten- 
tial. In the Jordan frame action, Eq. ([5]), matter fields 
are not explicitly coupled to the scalar field, hence, satis- 
fying the weak equivalence principle. Conformally trans- 
forming the metric g^ v and redefining the scalar field tp, 
the geometric part of the action and subsequently the 
Einstein field equations can be brought into their general 
rclativistic form. In this Einstein frame, matter fields are 
explicitly coupled to the new scalar field. If F = const, 
the two metrics are equal up to the constant factor F , in 
which case matter fields arc minimal-coupled. Models 
with non-constant F are consequently referred to as non- 
minimally-coupled scenarios. In this paper, we adopt the 
Jordan frame description. 

With the freedom to redefine the scalar field tp 
in Eq. ([5]), reducing the number of free functions of ip to 
two instead of three, we parametrize the non-minimally- 
coupled models using the Brans-Dicke (36j representation 



F 



Z 



_ u(tp) 



(3) 



For the minimally-coupled models, we set F(tp) = 
Z(tp) = 1. We refer to [37| for the more general descrip- 
tion of the background and perturbation equations using 
F(tp), U(tp), and Z[tp). All of the modified gravity and 
dark energy models discussed in this paper, except for 
DGP and possibly the phenomenological modifications 
of gravity (see §11 Cj) . can be formulated by an action in 
form of Eq. ©. 



1. Minimally- coupled examples 

We describe minimally-coupled scalar-tensor models 
by the modified Einstein-Hilbert action Eq. ([2]) with 
F(tp) = Z(tp) = I. Variation of the action with respect to 
the metric g^ v and the scalar field <p yields the Einstein 
field equation 

+df,pd v p - ^g^d a <pd a tp (4) 



and the conservation equation 

Dtp = dU/dtp. 



(5) 



Furthermore, energy-momentum conservation implies 
V„T£ = 0, where = 25S m /Sg^/^. 

Minimally-coupled scalar-tensor models do not mod- 
ify gravity since the Einstein equations for the Jordan 
frame metric are not modified, i.e., the Einstein tensor 
relates to the energy-momentum tensor as in GR and the 
scalar field can be interpreted as an effective fluid con- 
tribution to the matter and radiation components. The 



0.0 



-0.2 



FDSB 
SUGRA 



DOOM 
PNGB 



CPL-1 
CPL-2 




10 



FIG. 1: Dark energy equation of state associated with the 
quintessence models in £111 A II the freezing models moti- 
vated by the fermion condensate model with dynamical super- 
symmetry breaking (FDSB) and by supergravity (SUGRA), 
where the scalar field rolls along the potential at early times 
but then slows down at late times; the thawing models with 
collapse at doomsday (DOOM) and with a scalar-field po- 
tential motivated by the Pseudo-Nambu-Goldstone Boson 



(PNGB), where the scalar field is frozen at 



-1 at early 



times and begins to evolve at later times; and the phenomeno- 
logical quintessence models with constant and time-dependent 
dark energy equation of state. 



background field equation in a spatially-flat FLRW uni 
verse are 

H 2 



3H 2 = ^ P +~— p rz + U, 
-2HH' = K 2 (p + p) + H 2 p' 2 , 



(G) 
(7) 



where here and throughout the paper, primes denote 
derivatives with respect to In a. Furthermore, the energy- 
momentum conservation of the scalar field and matter 
gives 



H 2 p" + {H H' + 3H 2 )p' 



P 



dU 
dtp ' 
-3(p + p). 



(8) 
(9) 



The energy density and pressure associated with the 
scalar field also satisfy Eq. ©. By using the Friedmann 
equation and the energy-momentum conservation, we ob- 
tain the equation of state 



Pip 



(Hp , ) 2 -2U 
(Htp') 2 + 2U' 



(10) 



which remains in the regime w v G [—1, 1] for U > 0. 

In this paper, we study four different analytic poten- 
tials of tp and in addition, two potentials phenomenologi- 
cally defined by an effective Hubble expansion. We show 
the scalar-field equation of state for the different models 
in Fig. [U 
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a. Freezing models: As our first examples of 
quintessence potentials, we consider two freezing mod- 
els, in which the field rolls along the potential at early 
times but then slows down at late times (see, e.g., (HI), 



M 4+fi 
FDSB : [%) - 



<P 



n>0, (11) 



SUGRA : U(ip) 



M, 



4+n 



e o/ 



(12) 



where Mo, a, and n are the model-specific parameters. 
The potential in Eq. is motivated by the fermion con- 
densate model with a dynamical supcrsymmetry break- 
ing (FDSB) (H-lHI and the potential in Eq. dT2j) may 
arise from supergravity (SUGRA) models [42(. We set 
the FDSB potential with n = 6 in Eq. (fTTj) and the 
SUGRA potential with n = 11 and a = 1/2 in Eq. (fT2"]l. 
The scalar field ip at early epoch and the characteristic 
mass Mq are adjusted to match the present-day Hubble 
parameter H 2 (a = 1) = Hq. 

b. Thawing models: The next two potentials of in- 
terest here are classified as thawing models. In these sce- 
narios, the scalar field is frozen at early times with w v ~ 
— I and begins to evolve at late times (see, e.g., j38|). 



DOOM 
PNGB 



U 4 
M 4 < 



A I, 



'"ip 11 
[cap) 



n>0 (13) 
(14) 



The potential in Eq. ((13)) was studied in [43J as a model in 
which p rolls down the potential and collapses at dooms- 
day (DOOM). The potential in Eq. (Hi]) is motivated 
assuming the presence of an ultra-light Pseudo-Nambu- 
Goldstone Boson (PNGB) 0. We set the DOOM po- 
tential with n = 1 and Vq = 3Hq in Eq. (fT3")) and the 
PNGB potential with a = 0.445" 1 in Eq. flU}. The 
characteristic mass Mq and the initial conditions for ip 
are again adjusted to match H 2 (a = 1) = Hq. 

c. Phenomenological quintessence models: The last 
two quintessence models we analyze here are purely phe- 
nomenological and are based on the Chevallier-Polarski- 
Linder (CPL) [4^, |4(| parametrization of the dark energy 
equation of state w v = Wo + (1 — a)w a . In particular, we 
consider the two cases 



CPL - I 
CPL - 2 



w v = w , 



Wto = -a. 



(15) 
(16) 



Note that an equation of state of the form w v ~ —a 
can be made compatible with current cosmological obser- 
vations by introducing anisotropic stress (47j . Here, we 
shall, however, restrict to cases with no anisotropic stress 
(see, however, ^III E[) . which is a natural consequence of 
the minimal coupling with a constant F (see £|III B 21) . 
For the reconstruction of the scalar field, we treat its con- 
tribution to the expansion history as an effective fluid 

and ([10]) 



with equation of state w v and use Eqs. 
to determine U and ip' . We illustrate the dark energy 
equations of state of Eqs. (fT5"|) and ([TBI in Fig. [TJ where 
w = -0.8 in Eq. (|T5|). 



2. Non-minimally-coupled examples 

We describe non-minimally-coupled scalar-tensor mod- 
els in the modified Einstein-Hilbert action, Eq. @, in 
their Brans-Dicke representation, Eq. ©. Variation of 
the action with respect to the metric g^ v and the scalar 
field ip yields the modified Einstein and scalar-field equa- 
tions, 



R H" ~ 7j9nuRj = ^ 2 T^y + V M <9„(^ - g^Dcp - g^U 

+ ^ (d^tpdvtp - ^g t _ ll/ d a ipd a p^j , 

(17) 

2^Dp = -R-±(^) da(pd ^ + 2^. 
ip dip \ip J dip 



(18) 



We further require energy-momentum conservation, 
V^T£ = 0. 

Note that non-minimally-coupled scalar-tensor models 
modify gravity since the Einstein equations for the Jor- 
dan frame metric are not in their general relativistic form, 
i.e., the Einstein tensor does not relate to the energy- 
momentum tensor as in GR. The background modified 
Einstein equations and scalar-field equation are 



3ip H 



K 2 P +— if'' 

2 ip 



-3i?V + U, (19) 
-2ipHH' = K 2 (p + p) + H 2 -ip a 

+H 2 <p" 

+ (HH' -H 2 )p', (20) 

- \H 2 p" +(HH' + 'SH 2 )p'} =3(HH' + 2H 2 ) 

ip ' 

d f uj\ H 2 , 2 



dtp \p J 2 

dU 

dip 



(21) 



In addition, we have the usual energy-momentum con- 
servation of matter and radiation 



p' = -3(p+p). 



(22) 



As representative examples, we study two models of 
non-minimally-coupled scalar-tensor theories: f{R) grav- 
ity [481-52] or ui — Brans-Dicke models, i.e., with van- 
ishing kinetic term Z = in the Jordan frame action, 
Eq. ©, supplied with a scalar- field potential, and a 
scalar-tensor model with non-constant Brans-Dicke pa- 
rameter 1531] . 
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FIG. 2: Effective dark energy equation of state associated 
with different modified gravity models. As in Fig. [1] we 
consider six example models: two non-minimally-coupled 
scalar-tensor theories, f(R) gravity and the non-constant 
Brans-Dicke (NCBD) parameter model; the two branches 
of the DGP model, the self-accelerating (sDGP) and nor- 
mal (nDGP) branch; and two phenomenological modifications 
(PHEN-1 and PHEN-2). By design, PHEN-1 and the f(R) 
models considered here recover the background expansion his- 
tory of a ACDM universe, where w = — 1. Hence, their equa- 
tions of state overlap. 



a. f(R) gravity: The action of metric f(R) gravity 
models is obtained from Eqs. © and © by denning 



U = 
uj = 



= 1 + /b, 

Rf R -f 



o. 



(23) 
(24) 
(25) 



where f(R) is a nonlinear function of the Ricci scalar and 
fn = df(R)/dR. We focus here on designer f(R) mod- 
els [HI, where f(R) is reconstructed from a given Hubble 
expansion. We shall require a strict matter-dominated 
ACDM expansion history, i.e., we set H 2 — K 2 (p m +p^)/3 
in the f(R) Fricdmann equation, which is obtained by 
applying the replacements of Eqs. (|2"3")l . (|2"4"|) . and (|2~51 to 
Eq. (|19[) . The Fricdmann equation then becomes an in- 
homogeneous second-order differential equation for f(R) 
as a function of In a, 



1 



IT 



R' 



— T / = -H 2 (i-n m )— -. 

6H2 J 0V mJ H2 

This relation can be solved numerically with the initial 
conditions 



/(lnaO = AH 2 a\ - 6H 2 (l - Sl m ), 
f(\n ai )= P AH 2 a?, 



(27) 
(28) 

where p = ( — 7 + \/73)/4. A is an initial growing mode 
amplitude and characterizes a specific solution in the set 



of functions f(R) that recover the ACDM background. 
The amplitude of the decaying mode is set to zero in order 
to not violate high-curvature constraints [5J|- Instead 
of defining the solutions by setting A at an arbitrary 
redshift, we define them by a condition set today. This 
can be done by, e.g., either characterizing them at the 
background via fuo = fp(a = 1) or via the Compton 
wavelength parameter [54| 



B = 



1 + Ir h< 



(29) 



evaluated today, i.e., Bq = B(a = 1). GR is recovered in 
the case where Bq = or equivalently f^a = 0. 

Currently, constraints on f(R) gravity are of the or- 
der B < 1CT 3 or |/flo| (95% confidence level) 
from cosmological probes |55H57j . \fm\ < 10 -6 from So- 
lar System constraints [58| , corresponding to the require- 
ment of a chameleon-screened Milky Way halo [HI, [59[ 
(cf. [13]), and \fno\ 5, 10~ 7 from distance measurements 



and stellar physics [61(. An upper limit of \fm\ < 
is expected from future 21 cm survey data when com- 
bined with cosmic microwave background (CMB) data 
from Planck [62[. Also see 63 1 for expected constraints 
on f(R) gravity from a joint analysis of the bispectrum of 
the correlation of the integrated Sachs- Wolfe (ISW) effect 
with weak gravitational lensing and the CMB power spec- 
trum from Planck data. Despite the strong constraints on 
f(R) gravity, the model serves as a useful toy model for 
the study of the large-scale structure in modified gravity. 
We show the effective dark energy equations of state for 
the f(R) gravity models studied here in Fig. [51 The mod- 
els are designed to exactly match the ACDM expansion 
history, hence, w e f[ = — 1. 

b. Non- constant Brans-Dicke parameter model: Our 
second example with a non-minimal coup ling is equiva- 
lent to a scalar-tensor model studied in [53j. It has a 
non-constant Brans-Dicke (NCBD) parameter w(ip) and 
is defined by 



U 



w 



l-<p' 



(30) 
(31) 



where we set the two model parameters ujq = 1/40 and 
uq = —3 to emphasize near- horizon modified gravity ef- 
fects on the effective anisotropic stress in EIIIIB 31 For 
more details on the numerical integration, we refer to [53[ 
or Appendix [B] We illustrate the effective dark energy 
equation of state of the model in Fig. [5J which shows a 
phantom-crossing with w a e < — 1 and w e s > — 1 at early 
and late times, respectively. 



B. DGP braneworld gravity 

We extend our analysis of modified gravity to the well- 
studied DGP braneworld model 1 24j . Note that since the 
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PPF fits for the two DGP branches have been developed 
in H,[34|,[64|,[65|], we only provide a short schematic review 
and refrain from describing the modified Einstein equa- 
tions and their perturbations within DGP in the their 
full 5D form (see, however, (6i-[70jV 

In the DGP model, our universe is a (3 + l)-brane em- 
bedded in a 5D Minkowski space described by the action 



The DGP equation of state is given by 



S = --^ / d 5 x^g~R - JL / d^x^-gR 



2/i 2 



d A x 



(32) 



where 5D quantities are denoted by hats and 4D quan- 
tities are denoted by tildes. Matter fields, including a 
cosmological constant or brane tension, are represented 
by Lt and confined to the brane while only gravity can 
propagate through the full 5D bulk. We assume that 
there is no bulk tension. The constants k 2 and p 2 are 
proportional to the inverse Planck masses in the bulk and 
brane, respectively. Gravity on the brane is consequently 
modified at large scales. In particular, the crossover dis- 
tance r c = n 2 /2p 2 governs the transition from 5D to 
4D scalar-tensor gravity and determines the nonlinear 
screening, e.g., on scales smaller than the Vainshtein [2lJ 
radius r* = (j" 2 ^) 1 / 3 , nonlinear interactions return grav- 
ity to GR around a point mass with Schwarzschild radius 
V 

Variation of the action Eq. (I3"2l yields the modified 
Einstein equations on the brane, which for a homoge- 
neous and isotropic metric reduce to the modified Fricd- 
mann equation [66J 



H 



H 



3 ^' 



(33) 



where pi are the energy densities of various components 
on the brane and we assume vanishing spatial curvature 
on the brane. The sign a = ±1 defines the branch of 
the cosmological solutions. For a = +1, late-time ac- 
celeration occurs even without a cosmological constant 
A (6|| and so this branch is referred to as self-accelerating 
DGP (sDGP). For a = -1, DGP modifications slow the 
expansion rate and the branch is referred to as normal 
branch (nDGP). Here, a cosmological constant or addi- 
tional dark energy is required to achieve late-time accel- 
eration. For a spatially-flat matter-only universe with 
possible contribution of a brane tension or cosmological 
constant A, the Fricdmann equation becomes 



(34) 



with the density parameter associated with the crossover 
distance 



1 



1 — — £lf 



2H Q r c 



(35) 



where the second equality follows from the condition 
H(a = 1) = H . 



1 + IUDGP = 



M 2 Ei(! + w Qpi 

3H 2 + p 2 J2 iPi 



(36) 



where for cases with a cosmological constant, it is also 
useful to define the total effective dark energy p c ff — 
Pdgp + Pa and its equation of state 



1 + tUeff = (1 + WDGP) 



PDGP 



PDGP + PA 



(37) 



with pdgp = SH 2 / p 2 — EiPi- We show the effective 
dark energy equation of state associated with the DGP 
modifications in Fig. [2] where we use the fiducial cos- 
mological parameters and fl\ — or f2^ = 1.51 to de- 
termine the crossover distance in Eq. (j3"5")) for the sDGP 
or nDGP model, respectively. We use these parameter 
values throughout the paper. The effective dark energy 
equation of state in the nDGP scenario has a singular- 
ity at z ~ 1.5 where it jumps across the phantom divide 
with u> ff > — 1 and u> e ff < — 1 at early and late times, re- 
spectively. Note that the DGP background is completely 
defined by w e R and energy-momentum conservation of 
the brane components and the energy-density of the ef- 
fective fluid poff- 

The sDGP model is in conflict with cosmological data 
and can be ruled out at more than 5ct for both flat and 
nonflat universes [H|, [zj • If allowing for a brane tension 
in this branch (sDGP+A), a crossover scale of H$r c > 3 
(95% confidence level) is required for compatibility with 
cosmological observations, which essentially reduces it to 
a ACDM universe [H[ . A bound of the same order can 
be found for flat and nonflat nDGP universes [34[, where 
in the nonflat case, the cross correlation of the CMB 
temperature anisotropy with foreground galaxies through 
the ISW effect are used to break a degeneracy between 
the spatial curvature and the nDGP modification in the 
CMB temperature power spectrum and expansion his- 
tory [13, uM- Recently, constraints on a nDGP model 
with ACDM expansion history, i.e., with an appropri- 
ate additional dark energy contribution, have been ob- 
tained from redshift-space distortions [73[. Albeit the 
strong constraints and theoretical difficulties of the self- 
accelerating solution [74|, the DGP braneworld model 
still serves as a useful toy model for the study of the 
large-scale structure in modified gravity. 



C. Phenomenological modifications of gravity 

We can study deviations of GR in a more general 
context and allow for phenomenological modifications of 
gravity, where effective modifications in the Hubble ex- 
pansion as well as in the linear perturbations are intro- 
duced without a direct representation through a geomet- 
ric modification of the Einstcin-Hilbert action or Ein- 
stein equations but arc formulated by an effective fluid 
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contribution. Due to the energy-momentum conserva- 
tion of the matter component and the Bianchi identi- 
ties, the energy-momentum associated with this fluid is 
also covariantly conserved. While the linear perturba- 
tions of the fluid are characterized via the PPF formalism 
in £1111 E| at the background level, we describe the modifi- 
cation directly through its effective dark energy equation 
of state w c g(a) [cf. Eq. (j4"5"j) ]. The expansion history is 
then given by 



H 2 = y(Pm + Pcff), 



P'es = -3(1 



W e s)pef[, 



(38) 
(39) 



where we can write the dark energy density associated 
with the modification as 



PcS 



PcffO fl 

Peso exp 



3(l+t2e£f) 



1 l + ^cff(a') dQ / 



(40) 



Here, we implement the CPL paramctrization, which we 
have also used in £111 A 1 c[ i.e., 

w c e(a) = w + (1 - a)w a . (41) 

In this case, Eq. ([4T)|) simplifies to 

PcS = Pcsoa- 3{1+Wo+W ° } cxp[3w a (a - 1)]. (42) 

We shall consider two specific cases of phenomcnolog- 
ical modifications, 

PHEN - 1 : w = -1, w a = 0, (43) 
PHEN-2: w = -1.15, w a = 1.17. (44) 

The former follows a ACDM expansion history and the 
latter is motivated by the best-fit parameters to current 
cosmological observations found in [13| (see TOE]) . At 
the background level, the phenomenological modification 
of gravity is equivalent to a corresponding phenomeno- 
logical quintessence model and completely defined by 
the effective dark energy equation of state, Eq. (|4"Tj) . 
along with the Friedmann equation, Eq. (|38[) . and energy- 
momentum conservation, Eq. (f3T))) . In particular, PHEN- 
1 is equivalent to a ACDM model. They will only be 
completely characterized and distinguishable at the per- 
turbation level, which shall be elucidated in £|IIIEI Note 
that without a proper consideration of perturbations, the 
time-evolving dark energy models are gauge-dependent 
and theoretically inconsistent. 

We show the effective dark energy equations of state 
w e ft for our phenomenological modifications in Fig. [5] 
The first phenomenological modified gravity model ex- 
pands equivalently to ACDM, overlapping with the w of 
the designer f(R) model. The second phenomenologi- 
cal modified gravity model has a phantom crossing with 
w c s < — 1 at late times and approaches a matter-like 
equation of state at early times, i.e., w c s — 0. 



III. PARAMETRIZED POST-FRIEDMANN 
FRAMEWORK 

In 'analogy' to the parametrized post-Newtonian for- 
malism in the Solar System [l[, the PPF framework (see, 
e.g., [(j Il2| - [l4| ) provides a unified description of the lin- 
ear perturbation theory around the FLRW background 
for generalized modified gravity and dark energy mod- 
els. In this framework, the extra terms appearing in the 
Einstein field equations due to scalar-field contributions 
or modifications of gravity, e.g., Eqs. (@| and (JT7J) , are 
viewed as an effective dark energy fluid component, de- 
fined as 



1 cS — K u 1 m J 



(45) 



with energy-momentum conservation, V^T^g = 0, due 
to the Bianchi identities and energy-momentum conser- 
vation of the matter components. Thus, the usual cosmo- 
logical perturbation theory may be applied with each four 
degrees of freedom in both the perturbation of the metric 
and of the energy-momentum tensor. The Einstein and 
conservation equations fix four degrees of freedom and 
the gauge choice fixes another two. The remaining two 
degrees of freedom are then specified by two closure rela- 
tions, which are defined by the particular modified grav- 
ity or dark energy model. These two closure relations 
can be designed such that the effective fluid mimics the 
relations between the metric and matter perturbations 
given by the full perturbations of a particular modified 
gravity or dark energy cosmology. 

This simple and generalized treatment allows for an 
efficient and consistent computation of the evolution of 
perturbations on large scales in modified gravity and 
dark energy theories. Paramctrizations of this kind can, 
e.g., be applied to efficiently explore the parameter space 
of modified gravity and dark energy models employ- 
ing computationally intensive generalized Boltzmann lin- 
ear theory solvers and Markov chain Monte-Carlo tech- 
niques fl"ol - fl2l E3, [l! [zi|. A great advantage of the 
PPF formalism is that the corresponding numerical codes 
need only to be adapted once to employ the PPF modi- 
fications rather than for each non-standard model sepa- 
rately p^[75l - l77l |. Such a generalized approach motivates 
the application of the PPF formalism to our study, pro- 
viding the base upon which future horizon-scale tests of 
gravity and dark energy models employing, e.g., a multi- 
tracer analysis of galaxy-redshift survey data can be per- 
formed. 

Developing a PPF formalism for general theories of 
modified gr avity has been a subject to intensive study 
(see, e.g. [2T-fl8j|). In this paper, we follow @, adopt- 
ing a phenomenological approach to the PPF formalism 
(cf. [Ij]). It provides an efficient general framework to 
account for linear perturbations of modified gravity the- 
ories, consistently describing the sub-, near-, and super- 
horizon scales. Within this framework, paramctrizations 
of, e.g., f(R) and DGP gravity have been developed in Q 
and [a |34|, l&4|, respectively. In the case of DGP and 
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f(R) models (with large values of |/ro|)j the importance 
of the consistent inclusion of near-horizon modifications 
for the accurate determination of the low multipoles of 
the CMB temperature fluctuations have been pointed 
out in miMMIzlIzi ( cf - @)- Hence, in order to 
get consistent CMB constraints on the DGP and f(R) 
models, the DGP and f(R) PPF functions, respectively, 
have been employed in parameter analyses of these mod- 
els [lil, H(| [lH . The PPF formalism has also been utilized 
to consistently allow for phenomenological horizon-scale 
modifications of gravity and study constraints on such 
deviations from ACDM [47| . Thereby, the simplicity of 
the PPF formalism described in Q proved to be of great 
advantage in the efficient exploration of the parameter 
space through Boltzmann linear theory solvers [7l|, |zz| 
that have been employed to obtain accurate model con- 
straints. 

We briefly review the PPF framework of @ in mil Al 
In addition to the known PPF model descriptions in this 
formalism, in £1111 B 2\ we describe a simple procedure 
to calibrate near-horizon corrections to the sub-horizon 
PPF modifications for generalized minimally- and non- 
minimally-coupled scalar-tensor gravity models of the 
kind described by the modified Einstein-Hilbcrt action 
in Eq. ©, disclosing consistent near- horizon tests of 
these models through the CMB or multi-tracer analyses 
of galaxy clustering. We discuss the PPF fitting functions 
for f(R) gravity and the DGP model in flnCjand CDTDl 
respectively, and define the phenomenological modifica- 
tions of gravity in SjEIEj 



A. PPF formalism 

We introduce the linear PPF perturbation formalism 
of Q to compute linear perturbations of modified grav- 
ity and dark energy models. Thereby, we follow the lead 
of [1, H3, [U ■ For simplicity, we further restrict the lin- 
ear perturbations to a spatially-flat universe. The formal- 
ism can, however, easily be extended to include radiation 
components and spatial curvature [9| . In the longitudinal 
gauge for the metric perturbations (B = = Ht = 0, 
$ = Hi, *S? = A in Bardeen's notation [80j). assuming 
that there is no anisotropic stress n m = 0, the Einstein 
equations may be combined to yield @ 



matter 



4>_ 



2H 2 k 2 H 



PcffII c ff, 



2 ff2fc2 [PmA m + PcffAeff 

v eS -v m , 



+3(p c ff +Pcff)- 



Pcfrn cff 



(46) 



(47) 



where $+ = ($ + *) /2, $_ = ($- *)/2, k H = k/(a H) 
and w e s = Pcs/PcS is determined by the background 
relations. Using the energy-momentum conservation of 



A' m = -k H V m -3C, 

v' = k„y-v m , 



(48) 
(49) 



and the above Einstein equations, we obtain the velocity 
of the effective fluid 



2H 2 



K 2 (p c S +Poff) 



(50) 



and its pressure fluctuation 



Apcff = PeS A ofr - i/JeffAe ff - (p ff +Poff ) ^ + ('j , 

(51) 

where both quantities are in the longitudinal gauge and 
the comoving-gaugc curvature is £ = $ — V m /kn- These 
equations in fact define the effective fluid at the pertur- 
bation level. 

However, modified gravity theories are not character- 
ized by these properties of individual fluid components 
but are characterized by modifications of the Einstein 
equation. In the PPF formalism of @ , the Poisson equa- 
tion is modified as 



fc 2 ($_ 



(52) 



Together with the Einstein Eqs . (|46| ) and (|47|) . this yields 
the two PPF closure relations [6| 



2k 2 



T(a,k) = peff A ff + Pcffn e ff 



-3(p ff + Pes) -^i — — , (53) 



2H 2 k 2 

^g(a,fc)$_ = p cS U eS , 



where the metric ratio is defined as 



g(a,k) 



$ + $ + * 
$_ $ - * ' 



(54) 



(55) 



Under the PPF framework, the linear perturbations of 
modified gravity theories and dark energy models are 
completely described by the two additional functions 
r(a, k) and g(a, k), each of which defines the modification 
of the Poisson equation in Eq. (|53|) and the anisotropic 
pressure of the effective fluid in Eq. ([53]). respectively 
(cf. d-d !, EHH). In terms of these two PPF func- 
tions, the conservation equations become 



A' m = -k H V m -3(g + l)$'_ -3(l- 5 + </)$_ 

+3^ 
Vi = -V m + (g-l)k H *-, 



(56) 
(57) 



where we rewrite the metric perturbations in terms of 
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Given the metric ratio g(a, k) of each modified grav- 
ity model with its background evolution, the other PPF 
function T(a, k) is constrained in two limiting regimes, 
the super- and quasistatic sub-horizon scales. Since the 
energy-momentum is locally conserved even in modified 
gravity theories, the comoving gaug e curvature is con- 
served on super- horizon scales [Hj80j83l (cf. (54l.l68l.l83l]) 



C = o(k 2 H o, 



(58) 



and this relation constrains the PPF function T on hori- 
zon scales. With Eqs. (|5"2")) . (f5B|) . and ([ST)) , the relation 
on super-horizon scales in Eq. (|58[) can be used to yield 
the constraining equation on T in the limit kn — > as 



r 



5. 



(59) 



where the source function is 

1 H f 3 UqQh 



+ 



1 H 

-2.9 



2 H 2 cr 



(i + /c) 



k H 



.9+1 



(60) 



In the source function S, we parametrize the leading or- 
der correction to £' as 



lim -3- 



o(0- 



(61) 



For an effective scalar-tensor theory where photon 
geodesies are not affected by the additional scalar de- 
gree of freedom, the other constraint on the PPF func- 
tion r arises from the modified Poisson equation in the 
quasistatic regime on sub-horizon scales 



2(1 + fa) 



a 2 p„ 



X 



(62) 



which implies T — > fc&- in the limit ku — > 00, where the 
modification fc( a ) is defined at the linear perturbation 
level, ignoring nonlinear suppression mechanisms. Given 
the constraints in the two limiting regimes, a smooth 
interpolation of the PPF function rYa, k) can be achieved 
by solving the evolution equation Q 

(1 + c 2 k 2 H )[V + T + c 2 k 2 H {T - / G $_)] = S, (63) 

where the additional parameter cr relates the transition 
scale between the two limiting regimes to the Hubble 
scale. 

Therefore, the PPF function T(a, k) is further 
parametrized by the transition scale cr and two time- 
dependent functions /f(a) and /c(a), each of which re- 
lates the matter to the metric on super-horizon scales 
and defines the modified Poisson equation in the New- 
tonian regime, respectively. With the interpolation of 
Eq. (|6"5|) . the PPF formalism completely defines the lin- 
ear perturbations given the expansion history H or ui c ff , 
along with the PPF functions and parameter g, fe, fa, 



and cr- These functions can be described for modified 
gravity theories, e.g. for P GP @, [M HI (see fllHDj) . 
f(R) gravity [f| (see ijIIIC|) . and scalar-tensor models of 
the kind described by Eq. ^ (see ^III B 2|) . to obtain the 
linear fluctuations for the specific model. In order to de- 
termine the PPF linear perturbations, we need to solve 
the coupled system of differential and constraint equa- 
tions given by the conservation equations, Eqs. (|56[) and 
(|57|) . the evolution equation for T, Eq. (|63|) . with source 
Eq. and the modified Poisson equation, Eq. (|5^|) . 

See Appendix IA II for more details on the integration. 
Note that in mil B[ we reparametrize T in terms of a 
multiplicative correction to the standard Poisson equa- 
tion, i.e., S(a, k). In the limit of g = = fa = and 
correspondingly S = 1, the PPF perturbation equations 
recover the ACDM relations. 

Finally, the metric ratio g(a, k) can be parametrized 
by noting that g(a,kn = 0) = gsn(a) and g(a,kn = 
00) = <7qs(a) are scale- independent in the two limiting 
cases as m 



g(a,k) 



gsH(a) + (c g k H ) n o g QS (a) 
1 + {c g k H ) n * 



(64) 



where c g (a) and n g are the interpolation parameters that 
can be adjusted for each modified gravity theory. For 
more details, we refer to @ and Appendix IA 11 



B. PPF for scalar-tensor theory 

PPF fits for minimally- and non-minimally-coupled 
scalar-tensor theories or quintessence and general- 
ized/extended quintessence models, respectively, with 
horizon-scale consistency within the framework of @ 
have not been developed previously (see, however, [l4| for 
a description of scalar-tensor models in a different PPF 
approach). We shall therefore devote this subsection to 
introduce a calibration method for obtaining consistent 
near-horizon perturbations within general scalar-tensor 
models of the form of Eq. © . We then discuss the per- 
formance of this approach in comparison with the full 
model-specific fluctuations of the particular scalar-tensor 
theories. 

For our calibration method, we shall utilize the model- 
specific linear perturbation theory for scalar-tensor grav- 
ity and quintessence models, which we review in Ap- 
pendix [AJ Here, we discuss the corresponding quasistatic 
sub-horizon limit thereof, which we use to determine 
<7Qs(a) in Eq. ([M]) and the scale dependence of the metric 
ratio g(a, k). 



1. The quasistatic sub-horizon limit for linear perturbations 
in scalar-tensor gravity 

In the quasistatic approximation at sub-horizon scales, 
the perturbed modified Einstein equations of non- 
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minimally-coupled scalar-tensor models in their Brans- 
Dicke representation (see Appendix IA 31) can be simpli- 
fied to obtain the relations (see, e.g.. [37ll83|) 



9 



1 



2uj + 3 a 2 M 2 + k 2 ip ' 



fc 2 $_ 



—a p m A m . 

2ip 



(65) 
(66) 



In the limit k ^> aM^ -1 / 2 , Eq. (|65|) becomes scale- 
independent with the metric ratio g — > — l/(2w + 3). The 
mass of the scalar field M, governing the scale depen- 
dence in Eq. (|65j) . is determined from Eq. (fl"8|) . replacing 
R with the trace of Eq. (fT7|). i.e., 

(2w + 3)Qp = K 2 T - cj v (a Q ^) 2 - 4U" + 2<p U v . (67) 

By defining the effective potential Dip = V e s,ip, the mass 
of the scalar field is 



model-specific linear perturbation theory defined in Ap- 
pendix IA 31 for each non-minimally-coupled scalar-tensor 
model at a given fci, which can be chosen through the 
required range of applicability of the fit. In this pa- 
per, we shall work with k\ = Ho- The resulting met- 
ric perturbations can be used to obtain the metric ratio 
5i(o) = g(fl,k\), which then yields 



5sh(o) = ffi(o) 



9i{a) 



1 



2uj 



9 



VaM 



(74) 



In order to relate the metric to the matter perturba- 
tions, we reparametrize the modification to the Poisson 
equation of the lensing potential in Eq. (|52l) by [H 



fc 2 $_ = E(a, fc)ya 2 p m A m , 



where we have defined 



(75) 



M 



dip' 



off • 



(68) 



Together with the conservation equations, Eq. (|A14I) . 
which becomes = -k H V m ((' = 0), and Eq. (|A15D . 
Eqs. (jnSJ) an d with initial condition A mj ; fully define 
the quasistatic sub- horizon linear perturbations, i.e., 



a;; +(2 

,2 i\j2 i ( o. , i /|\;„2,„ tt2 



— ) A' 



3 (2w + 3)a 2 M 2 + (2cj + 4)k 2 ip H$Q, n 
2^ (2w + 3)(a 2 M 2 + k 2 ip) H 2 a 3 



A m = 0. 



(69) 



2. PPF horizon calibration 



Based on the quasistatic sub-horizon non-minimally- 
coupled scalar-tensor perturbation limit for the met- 
ric ratio, Eq. (|65|) . and the PPF interpolation formula, 
Eq. (|64|). we extrapolate a PPF fitting function for g 
for general scalar-tensor theories of the kind considered 
in Eq. ([3J in Brans-Dickc representation. The quasistatic 
sub-horizon part of our scalar-tensor PPF fit for the met- 
ric ratio is obtained from Eq. (|65[) , hence, we define 



1 



9QS 



2w + 3 

c 9 = SP W 
n g = 2, 



(70) 

(71) 
(72) 



in Eq. which yields the scalar-tensor PPF fit for the 
metric ratio 

, M 1 a 2 M 2 (2u; + 3)g sli (a)~ipk 2 
9i a > k ) = o 7^ 277^ 72 > ( 73 ) 



2a; + 3 



a 2 M 2 + p k 2 



where the mass of the scalar field M is obtained from 
Eq. (|fJT)l . We determine <7sH(a) by solving the full 



E(a, k) = 



1 + r/$_' 



which reduces to 



E(a, k) 



= E 



G 



QS 



(76) 



(77) 



in the quasistatic sub-horizon limit, recovering Eq. ([62 
This sets the PPF parameter 



/g = <p ~ 1 



(78) 



from Eq. (|66[) for non-minimally-coupled scalar-tensor 
models, i.e., Eqs = y -1 . In order to obtain E(a, k) on 
super- and near-horizon scales, we use the full model- 
specific perturbations at k\ utilized in the determination 
of ffsn(a) in Eq. (|63|) and (|60)) to obtain f^(a). Note that 
we refrain here from using simple numbers or scalings for 
/f(a) based on c/sH~(a), providing a simple and general 
approach for determining /(-(a). 

The PPF function determined by this procedure 
together with fa, g, and w e g, which is obtained from 
the background equations, Eqs. (|T§| to (|2"2l . fully de- 
termine the linear PPF perturbations corresponding to 
the model-specific perturbations of the non-minimally- 
coupled scalar-tensor theory of Appendix IA 31 up to the 
transition scale cr, which remains the only free PPF pa- 
rameter that has to be calibrated to obtain the correct 
scale dependence in the modification of the Poisson equa- 
tion of the lensing potential, Eq. (|75|) . 

For minimally-coupled models during matter domina- 
tion, the metric ratio is 



9 = 0, 



(79) 



which follows from F — const and the Einstein equa- 
tion $+ = -5F/F = (see, e.g., |2Z|)- From this and 
Eq. (lA"5"]l. it follows that 



fa = 



(80) 



11 



in minimally-coupled models. However, in general, E ^ 1 
at near-horizon scales since dark energy may cluster at 
the largest scales. Note that for minimally-coupled mod- 
els, the Poisson equation should not be regarded as mod- 
ified in the sense of changing gravity but contributions 
of scalar-field perturbations are parametrized as an effec- 
tive modification to the relation between metric and mat- 
ter perturbations. By calibrating cr, we can apply the 
same procedure as described here for the non-minimally- 
coupled scalar-tensor theories to obtain /f(a) at k\ for 
the quintessence models and hence determine the linear 
PPF perturbations corresponding to the full quintessence 
perturbations of Appendix IA 21 



3. Performance of PPF fits 

We test our PPF fits for the metric ratio g(a, k) and 
the modification of the Poisson equation E(a, k) de- 
rived in ^IH B 2 1 against their predictions from the model- 
specific linear perturbations of the quintessence models 
described in Appendix IA 21 (see Fig. [3]) and the non- 
minimally-coupled scalar-tensor theories of Appendix IA 31 
(see Fig. [4]). For the cosmological and model parame- 
ters, we adopt the values defined in £JTX| In the case of 
f(R) gravity, we also give the original PPF fit of [HH 
(see ^III C|) in the first row of Fig. 21 which is in good 
agreement with the results from mil B 21 Note that due 
to our calibration of at fcj = Hq, our fit captures the 
'true' Efa, k) slightly better than the original f(R) PPF 
fits of Q, which use = — 5sh/3. We choose a rather 
large value of |/ro| to highlight differences between the 
different PPF predictions as well as the full perturbation 
theory for f(R) gravity. The deviations become smaller 
with smaller values of |/ao|- Note, however, that an up- 
per bound of I/roI < 0.35 and \f R0 \ < 0.07 (95% con- 
fidence level) is inferred from CMB only and in com- 
bination with galaxy-ISW data [56j|, respectively. For 
I /ho I = 0.05 chosen in Fig. [4] at the horizon today, the 
metric ratio can be of order g ~ 0.1, which is a factor 
of ~ 2 larger than the result we would obtain if ignoring 
the super- horizon contribution, i.e., setting g$n = in 
Eq. ([MD- 

Our PPF fits arc in good agreement with the 
model-specific perturbations of the minimally- and non- 
minimally-coupled scalar-tensor theories. Note that 
there remains potential to improve the correspondence 
in E(a, k) by introducing time- and/or scale dependence 
in the transition scale cr- 



C. PPF for f(R) gravity 

The PPF expressions for f(R) gravity were developed 
in Q and have been shown to provide excellent fits to 
the full perturbation evolutions while at the same time 
making the integration computationally less expensive 0, 
|56| . On super-horizon scales, i.e., in the limit of k R — > 0, 



the metric ratio gsn = $+/$_ is defined through the 
metric perturbations from solving 



1 



El 



B' 



1 - B 

El _ Ell 
IT ~ ~w~ 



B H ' 

B' 
1 - B 



$ = 0, (81) 



and using the relation 



1 B ' 



(82) 



where B is the Compton wavelength parameter defined 
in Eq. (|29|) . Eq. (|82|) follows from conservation of curva- 
ture fluctuation (£' = 0) and momentum, considering the 
super-horizon anisotropy relation $ + = —BH'aV m /k. 

In the quasistatic regime, gQs = — 1/3 and at interme- 
diate scales, g(a, k) can be determined through the in- 
terpolation formula Eq. (fM]) . where c g = 0.71y/B and 
n g = 2 | Further, f c = -£sh/3, /g = /fl and 
cr = 1 [6|. Note that the approximation breaks down 
at intermediate scales and low redshifts for values of 
^o^l [5(| . which are, however, strongly ruled out by 
observations I55l-l57l. 



Correspondence to scalar-tensor PPF 

f(R) gravity is equivalent to a non- minimally-coupled 
scalar-tensor theory with the specifications given in 
Eqs. (HU), (HH, and ([25]). Hence, the PPF fits for gen- 
eral non-minimally-coupled scalar-tensor theories con- 
structed in §111 Bl should reproduce the f(R) PPF de- 
scription in the corresponding limit. 

In /(i?) gravity, the scalar- field equation is obtained 
from setting lo = 0, <p = 1 + fn-, and U = (R/r — /)/2 
in Eq. (|67|). hence, 



af R =l[(l-f R )R + 2f + ^T]^^ 



(83) 



which defines fc = f R and the mass of the f R field 



At high curvatures, 



Mf 



1 1 

3~ 



Jrr 



f R 1 H R' 



3 B H' 



RR 



(84) 



(85) 



where for a background evolution recovering a ACDM 
expansion history, this becomes M"j R = 2H 2 /B with B 
determined by Eq. (|29p . For the interpolation weights, 
this yields 



(1 + fnh 



0.71 VB 



(86) 
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FIG. 3: Modification E of the Poisson equation, Eq. ((75}, relating the lensing potential to the matter density fluctuation, due 
to contributions of the scalar field (E — 1 in ACDM). We show the clustering of the scalar field for the six minimally-coupled 
scalar-tensor (or quintessence) models defined in AH A II at large scales. The quintessence field does not introduce anisotropic 
stress (g — 0). Given the background evolution, the two PPF functions completely define the quintessence perturbations. The 
PPF fits developed in TOI B 21 (dashed lines) provide good fits to the model-specific quintessence perturbations (solid lines). 
The transition from the quasistatic sub-horizon to the super-horizon limits, governed by the PPF parameter cr, is calibrated 
at ki = Ho, where the fit overlaps with the predictions from the full perturbations. Note that E does not recover the ACDM 
limit at early times for the CPL-2 model due to the strong background contribution of the scalar field at small a. 



13 



From left: H l k = 100, 50, 20, 10, 5, 2, 1 

|/ RO |=5xlO " 

- - - - PPF (ori 
PPF(c r = 




1.05 



1.04 



1.03 



1.02 



1.01 



1.00 



|/ R o|=5xlO- 2 
PPF (original) 
PPF(c r =l) 




In a 



-3.0 -2.5 -2.0 -1.5 
In a 



-1.0 -0.5 0.0 



NCBD 

PPF (c r =0.65) 




1.04 
1.02 
1.00 
3 0.98 
0.96 
0.94 
0.92 



NCBD 

PPF (c r =0.65) 



-3.0 



-2.5 



In a 




FIG. 4: Metric ratio g describing the effective anisotropic stress (left-hand column) and the modification of the Poisson 
equation E (right-hand column) for the designer f(R) gravity f i]II A 2 a|) (upper row) and the non-constant Brans-Dicke (NCBD) 
parameter model defined in Eqs. (|30p and (|3ip f ^II A 2 b|l (lower row). The calibration technique developed in ijHI B 21 produces 
a slightly improved fit for E compared to the original f(R) PPF fit and is governed by the PPF parameter cr- Note that for 
f(R) gravity with small values of \fno\, the degree of correspondence between the calibrated and original PPF as well as the 
'true' functions g and E improve. 



for \ffj\ <C 1 with n g = 2. Since uj = 0, we have gos = 
-1/3 in Eq. §5§. This recovers the PPF fits of at 
small and intermediate scales. 

We make comparisons of the super-, near-, and sub- 
horizon PPF fits in mil B 31 Note that at early times, 
we can neglect the time derivative in Eq. (|82[) . which 
leads to gsa — > —B/2. Using this approximation, we 
may obtain a simplified super-horizon PPF fit for g with- 
out solving Eq. (jSTj) through g SH ^ — (B / 2)Q m (a) n with 
n m (o) = {H /H) 2 n m a- 3 , where n £ (0.5,0.6) is found 
here to give good agreement with the full perturba- 
tions. Here, we shall, however, not use this approach 
but wc remind the reader that simple approximations 



of this form can be very useful in parameter estimation 
analyses emplo ying, e .g., PPF Boltzmann linear theory 
solvers [HI, [fj I75l477j in a Markov chain Monte- Carlo 
exploration of the cosmological parameter space to ob- 
tain constraints on the gravitational models. For this 
purpose, the perturbed field equations need to be solved 
several thousand times. The PPF approach can signifi- 
cantly increase the efficiency of the integration. For an 
application of the PPF framework in f(R) gravity model 
constraints, see, e.g., [Hfjj]. 
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D. PPF for DGP 

The PPF functions for sDGP were derived in Q (also 
see [HI) and have been shown to provide excellent fits 
to the full perturbation evolutions while at the same 
time making the integration computationally less expen- 
sive [1, l34l. FzTl] - They encapsulate bulk effects in an 
effective 3 + 1 description. In [34j ]. these results were 
extrapolated to cases with nonvanishing brane tension 
and curvature. The PPF functions for nDGP were de- 
scribed in [64[ based on bulk calculations of [69| and [70[ 
with brane tension but no curvature. Simplified fits with 
additional extrapolation for spatial curvature are given 
in 34[. Here, we adopt the fits of 34| assuming spatial 
flatness K = 0. 

In the quasistatic high-fc limit, the DGP model pre- 
dicts la. 16511 



9qs 



H' 

1 - 2ct H r c ( H 

3H 



(87) 



On super-horizon scales, the sDGP and nDGP fit is [H,|34| 



in Eq. ([64]) such that the modifications of the model are 
defined by 



g(a, k) = 



9oa 



k 2 H 



(90) 



where we set go = 10 for illustration. Note that in such 
a model, departures of ACDM disappear on sub-horizon 
scales, i.e., g < 10~ 3 at k > 10 2 Hq today, whereas g = 5 
at k = Hq. 

Our second phenomenological model (PHEN-2) is mo- 
tivated by the results of 47J • The expansion history is de- 
fined through the effective dark energy equation of state 
in CPL form, i.e., wo and w a in Eq. (|4"2"|) . We further de- 
fine /g = /c = and cr = 1. The metric ratio is chosen 
to be a constant, 



g(a,k) =g , 



(91) 



S§g GP («) 



.9sg GP («) 



9 



8H r c - 1 
1 



0.51 



H r c - 1.08 



2Hr r 



(89) 



respectively. For the interpolation on intermediate scales 
through Eqs. (|63|) and (|64[) . we have c g = 0.14 and cr = 1 
in sDGP and c g = 0.4 and c r = 0.15 in nDGP 0. 
Furthermore, n g = 3 and the PPF functions that relate 
the metric to the density can be set as /f = 0.4gsH and 

/ G = o| 

Note that since the PPF fits to the full 5D DGP per- 
turbations have been studied in [(| |64| and we do not 
repeat these calculations here, we refrain from showing 
any DGP predictions for g(a, k) and S(a, k) and refer 
to @, Hi] for comparisons of the PPF with the full DGP 
perturbations. 



E. PPF as phenomenology 

In order to study more general effects on galaxy clus- 
tering caused by deviations from the concordance model, 
in addition to the scalar-tensor and DGP models, we 
further construct two phenomenological modifications of 
gravity (or non-standard cosmological models) based on 
the PPF formalism described in fJTTTJ We assume for both 
models a standard quasistatic sub-horizon Poisson equa- 
tion or negligible modifications thereof, i.e., fc = 0, and, 
furthermore, /f = with cr = 1. 

In the first model (PHEN-1), we assume a ACDM ex- 
pansion history, i.e., wo = — 1 and w a = in Eq. (|42[) . and 
an effective super-horizon anisotropic stress described by 
gsn(o.) = goa 3 such that modifications vanish at early 
times. We further choose gqs — 0, c g — 1, and n g = 2 



i.e., modifications are also present at the largest scales 
and early times. We motivate our choice of wq, w a , and go 
by the best-fit parameters found in (47| . i.e., wq ~ —1.15, 
w a ~ 1.17, and go ~ —0.15. Remarkably, such a model 
was found to be consistent with current cosmological ob- 
servations on linear scales, simultaneously including all 
measured multipolcs of the CMB data, their cross corre- 
lation with foreground galaxies through the ISW effect, 
as well as geometric, weak lensing, and clustering probes 
while slightly increasing the best-fit likelihood and en- 
hancin g th e average likelihood over the posterior distri- 
bution [471 ]. The corresponding marginalized likelihoods 
were found to be consistent with ACDM [47j . The key 
feature of this strong non-standard consistency is a de- 
generacy between the early anisotropic stress and w a 
when w a tn wq in the CMB data [47[. Note, however, 
that other cosmological parameters like the matter den- 
sity O m arc not set to the best-fit values of |47| but are 
fixed to the same values that are used here for all models 
(see ^TT| and, hence, deviations from ACDM predictions 
may appear stronger. 

We stress that the phenomenological models given here 
may not be physical, i.e., it is not clear whether more 
rigorous theoretical models with these properties may 
be constructed. The PPF formalism, however, offers 
the possibility of testing such non-standard cosmologies 
while retaining energy-momentum conservation of the 
fluid contributions and avoiding gauge artifacts as ap- 
parent deviations from ACDM. Albeit the possibility of 
these phenomenological modifications to be unphysical, 
if ACDM is the correct model, we ultimately aim at rul- 
ing out such deviations through observations. Moreover, 
such modifications offer a great opportunity to test the 
constraining power of future surveys. 
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IV. SIGNATURES OF MODIFIED GRAVITY 
AND DARK ENERGY MODELS IN GALAXY 
CLUSTERING 

With the PPF fits for modified gravity and dark en- 
ergy models developed in mill we consistently compute 
the galaxy power spectrum from the sub-horizon to the 
near-horizon scales. In ^IV At we briefly review the rela- 
tivistic formula for the observed galaxy fluctuation field 
and discuss its applications to modified gravity and dark 
energy models. In §IVB1 we estimate the detection sig- 
nificance of those models expected with future galaxy 
surveys. 



A. Relativistic effects in galaxy clustering 

The relativistic description of galaxy clustering has 
been developed in the past few years [25M3ll [33| . Given 
the observed redshift and the observed galaxy position 
on the sky, the full relativistic formula for the observed 
galaxy fluctuation can be derived by tracing the pho- 
ton path backward in time. Since the observed quan- 
tities and the photon path are affected by not only the 
matter fluctuations but also the relativistic contributions 
of the gravitational potential, the relations among the 
source, the observer, and the metric perturbations are 
nontrivial and they require a relativistic treatment for 
solving the geodesic equation. In this manner, the full 
relativistic description provides a complete picture of all 
the effects in galaxy clustering such as the redshift-spacc 
distortion, the gravitational lensing, the Sachs- Wolfe ef- 
fect, and their relativistic corrections. Furthermore, this 
procedure is necessary because all these effects are not 
individually separable and what we measure in galaxy 
surveys are the number of galaxies, their observed red- 
shifts, and the positions on the sky, not the individual 
effects. We refer to [2^, [2(| for the derivation and subtle 
gauge issues associated with it. 

The matter and metric perturbations along the pho- 
ton path affect the observed galaxy fluctuation. Fur- 
ther contributions arise from the vector and the tensor 
perturbations. In generic modified gravity theories, the 
vector perturbations decay with time as in the standard 
ACDM model. Furthermore, the primordial tensor per- 
turbations are well constrained by the CMB observations. 
Therefore, in modified gravity theories and dark energy 
models that mimic the ACDM model at early times, we 
can safely ignore the vector and the tensor contributions. 
With only the scalar contributions, the observed galaxy 
fluctuation can be obtained as [III |33| 

A° bs = bA m -eSz v + y + 23> + V + 35z x 



-2K-H- 



6z x 
dz \ a H 



5 P 8V L ,(92) 



where b is the linear bias factor, the evolution bias factor 
e = d In fig/d ln(l + z) is the redshift evolution of the ob- 
served galaxy population with the galaxy number density 



n g , p = — 0.4dlogn 9 /dlogL is the slope of the luminosity 
function, Sz v and Sz x are the lapse in the observed red- 
shift in the matter comoving and the longitudinal gauges, 
V is the line-of-sight velocity, 5r x and K, are radial and 
angular displacements, and 8T>l is the fluctuation in the 
luminosity distance. Hereby, the subscript v and x de- 
notes the dark matter comoving (v = 0) and conformal 
Newtonian (zero-shear x = 0) gauge, respectively. 

In addition to the modulation of the galaxy number 
density due to the matter fluctuation &A m , the physical 
origin of the numerous terms in Eq. (|92[) are the vol- 
ume and source distortions. Since we define the observed 
galaxy samples in terms of the observed redshift and flux, 
the fluctuations in the observed redshift and the flux give 
rise to the source effects eSz v and 5p5T>L- The metric 
potential and velocity terms define the Lorentz frame of 
the source galaxies. The remaining terms in Eq. (|92[) 
arise due to the volume distortion, as the volume element 
we assign in observations based on the observed redshift 
and angle is different from the volume occupied by the 
observed galaxies. The volume distortion is decomposed 
as the radial Sr x and the angular K, displacements and 
the distortion in the observed redshift. The notation is 
written in a physically straightforward manner. 

The application of the relativistic description to the 
modified gravity theories and dark energy models is 
straightforward. Though the full relativistic formula in 
Eq. (|92p is obtained by solving the photon geodesic equa- 
tion in GR, the sole assumption that leads to Eq. (f9"2")l 
is that the space-time is described by the FLRW metric 
and photons follow the geodesic. This assumption is valid 
for all modified gravity and dark energy models consid- 
ered in this paper. Therefore, Eq. (j9"2")l is applicable to 
those models without modification. We further simplify 
Eq. (f9"2"| for the power spectrum analysis by ignoring the 
complication of survey geometry and projected contri- 
butions to A° bs . The former can be readily dealt with 
once the survey boundaries are specified while the lat- 
ter is negligible in the power spectrum analysis (2|| [Hf . 
In our metric notation, the observed galaxy fluctuation 
becomes 



. obs 



b A m - n\k H V m - - - {hp - 2) $ 



e + 1 + ^ + (5p - 2) (l - 

H' f 1 

e + — + (5p - 2) I 1 - 



1 



aH r 



H 



a H r 



(93) 



where /x& is the cosine of the angle between the Fourier 
mode and the line-of-sight direction. Note that we do 
not equate the two metric potentials "J and — <& as they 
differ in modified gravity scenarios (g 7^ 0). 

The galaxy fluctuation in Eq. (|93p is written in terms 
of metric perturbations in the longitudinal gauge. Note 
that since the full relativistic formula in Eq. (|9"2"|) is gauge- 
invariant, it can be evaluated at any choice of gauge con- 



16 



0.48 : z=0 
0.47 1 

0.46; 

\ 

ft. 0.45 1 7 ^ 
0.44 1 
0.43 ] 
0.42 1 



FDSB 1, 0.55 
SUGRA 
DOOM 

PNGB j o.50 
CPL-1 



0.45 



0.40 



CPL-2 
ACDM 



10 

k/H Q 



20 



50 



100 




NCBD 
I/roI=10" 
I/roI=10- 
I/roI=10- 



10 



\fw>\=W~* 

- I/roI=io- 5 

- I/roI=1(T 6 
ACDM 



50 100 



500 1000 



k/H 
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the Newtonian growth rate T reduces to the usual logarithmic growth rate / of structure in Eq. (|98p . The gray-shaded region 
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ditions, including the longitudinal gauge in Eq. (|93|) . The 
relation of the metric potentials to the matter density 
fluctuation A m and its velocity V m depends on the par- 
ticular gravitational model. 

We shall first consider the case of a ACDM universe. 
Using the Einstein equations and the conservation equa- 
tions in GR, the galaxy fluctuation in Eq. can be 



scaled with the matter density fluctuation as [26J, |31|, |3c 
A° bs = (94) 



a Newt. , J^a _ ■ ^ a 
T f.2 ^ m ^ k kil 



where the two additional terms 
3 



2 

+ef 



n 



H' ( 
+- — + (5p - 2) I 1 - 



H 



aHr 



(95) 
(96) 



represent the gravitational potential and the velocity 
contributions to A° bs . Here, the scale dependence of 
these contributions is explicitly removed and we have as- 
sumed a ACDM universe. Moreover, we have defined 



i /H 2 . The galaxy clustering in red- 



shift space is often described by the Newtonian contribu- 
tion l8al 



.Newt, 
g 



bA m - f ilk H V m = bA m + f f i 2 k A m , (97) 



where the logarithmic growth rate of structure is 



f = 



d In A m 
din a 



A' 
A 



=^c±-k H ^. (98) 



The last relation is derived from the conservation equa- 
tion in Eq. (|48|) and becomes exact in a ACDM universe 
with GR, where (' — 0. It is shown [33[ that with 
the hindsight of the full relativistic formula, the veloc- 
ity terms in Eq. (|92[) can be derived by using Newtonian 
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dynamics but we refer Eq. ([97)1 to the Newtonian contri- 
bution because it is the only non-negligible contribution 
in the Newtonian limit kn — > oo. 

In a ACDM universe with GR, the key ingredients /, 
V , and TZ in Eq. ([94]) are all scale- independent, simply 
because the scale dependence of the relation among VP, 
$, V m , A m is explicitly taken into account and removed. 
However, these coefficients in Eq. (|M|) become scale- 
dependent in modified gravity and dark energy models. 
To properly account for these modifications, we first de- 
fine the Newtonian growth rate as the ratio of the velocity 
to the matter density 



-k 



Kr, 



H- 



(99) 



and the conservation equation in Eq. (|48|) yields the log- 
arithmic growth rate of structure 



c 



T. 



(100) 



In the Newtonian limit, where the effective fluid is well 
approximated as a single matter fluid with vanishing 
pressure, the matter comoving-gaugc curvature is con- 
served, Q' — 0, and both quantities are identical, / = J-, 
providing a further justification of its terminology. How- 
ever, it is really the ratio of the velocity to the matter 
density and the redshift-space distortion in Eq. (|97|) con- 
strains J 7 , not /. The logarithmic growth rate / can be 
measured, for example, in the cluster abundance analysis. 
However, for the models we consider in CTfl the comoving- 
gaugc curvature can vary in time on large scales and, 
hence, the two growth rates / and T are different. The 
deviation due to the comoving curvature is described by 
the effective fluid as m 



2H 2 



a 2 (p eS +Pes) 



(101) 



By using Eqs. {52$, {TSJ, and dMJ), we derive the rela- 
tivistic description of galaxy clustering in a PPF universe 
based on Eq. (jM)) . 



T 



2k%f + 9n m (a) [g 1 + (g + !)(£'/£ + /) - 2g] £ 



2k 2 H - 6H'/H 



V = eT-^n m (a)lg'-2g+(l+g) f^ + f 



+(5p-2)(l + g) + (l-g) 



e-f4 + (5p-2)ll 




(102) 

I 

(103) 
(104) 



In Figs. [S] and [51 we show the predictions for T for the 
different dark energy and gravitational models. Fig. [5] 
illustrates the corresponding model predictions for the 



rclativistic corrections, i.e., the velocity and potential 
coupling, TZ and V, respectively. We adopt the model 
and cosmological parameter values used in £jll] and the 
perturbations are obtained employing the linear PPF 
framework discussed in mill Note that at fixed redshift, 
TZ cx J 7 , where for e = 3 and p = 0.4, we have 

3 

TZ = - [2 - fl m (a) - (1 + w cS + w' cS \na)Q cS (a)} T 

(105) 

with f2 ff(a) = H 2 K 2 p e ff/3. Given this proportional- 
ity, simply being the consequence of TZ describing the 
coupling to the velocity knV m in Eq. (|9"3"]l and JA„, = 
— kjjV m , we can directly relate the signatures and cor- 
responding measurement significances for the modified 
gravity and dark energy cosmologies in J- and TZ. 



B. Measurement significance 

In a multi-tracer analysis (87j of galaxy-redshift survey 
data, when the observed galaxies are divided into differ- 
ent samples of different values of galaxy bias the obser- 
vation of modes transverse to the line-of-sight direction 
(pk = 0) yields a measurement of the relation between 
the different galaxy biases, i.e., b t up to the normaliza- 
tion b, as is evident in Eq. (|9"T|) . Neglecting stochastic- 
ity between the dark matter and the galaxy fluctuation 
fields, subsequent observations of a mode at ^ give 
a measurement of which is free of sample variance. 
One can furthermore determine the bias, free of cosmic 
variance, by cross-correlating galaxy-redshift survey data 
with weak gravitational shear fields [H| . A combination 
of these two methods therefore, in principle, yields a mea- 
surement of the logarithmic growth rate / which is free 
of sample variance [89[ . 

In practice, however, galaxies are discrete objects, 
and their measurements necessarily involve intrinsic shot 
noise, let alone other errors associated with the measure- 
ments. Once the sample variance limit is eliminated, the 
shot-noise errors are the dominant uncertainties. Since 
the dark matter halos in which the galaxies reside have 
finite size in their extent, their sampling errors deviate 
from the Poisson noise of point sources. Based on this 
observation, proper mass-dependent weights can be ob- 
tained to approximate dark matter halos as the dark 
matter distribution, reducing the stochasticity between 
them. This shot-noise canceling technique [90l - fQ3l | can 
be combined with the multi-tracer method to maximize 
the advantage of both methods. 

Bernstein and Cai [94| performed a Fisher-matrix anal- 
ysis to derive the cosmological constraining power from 
a large-scale galaxy survey combined with weak lensing 
measurements. Adopting the Kaiser formula in Eq. (|9"T)) 
and applying the multi-tracer method with the shot- 
noise canceling technique, they find that in a future 
half-sky galaxy-redshift survey, including halos of mass 
M > lO^h' 1 M Q and modes with k < 0.03 h/Mpc at 
z = 0.5 ±0.1, assuming that weak lensing can determine 
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0.0 0.0 
1.6 0.1 



TABLE I: Example models of gravitational and quintessence theories studied in ijll]with their corresponding PPF functions 
and parameter describing the linear perturbations (see ^jlllj) . Hereby, 'cal.' refers to the application of our calibration technique 
l ^III B 2[) . In the last column, we present approximate significance levels of the deviations of the modified gravity and dark 
energy models based on the constraints that may be obtained in a multi-tracer analysis of future galaxy survey data at 
z — ~ 1 (see ^IV B[) . For simplicity assuming fixed cosmological parameters and scale-independence in J 7 , 1Z, and V, the 
measurement significance of the relativistic correction TZ is comparable to its counterpart for T, which in ACDM corresponds to 
the logarithmic matter-density growth rate /. Thus, relativistic corrections to the Newtonian description of galaxy clustering 
need to be taken into account for deriving consistent model constraints in a multi-tracer analysis of future galaxy-redshift 
survey data. This conclusion holds even in the case of including more small-scale modes in the analysis, which increases the 
measurement significance of deviations in T by a factor of ~ 6. Note that for f(R) gravity, we have used |/ijo| = 10 _1 . Other 
model parameters are given in fjll] 



the bias to an uncertainty of <Jb = 0.01, the growth index 
7 can be measured to an uncertainty of about 9% 94j . 
where / ~ f2 m (a) 7 . Note that if including smaller scales, 
i.e., A; < 0.1 /i/Mpc, assuming that the linear Kaiser for- 
mula still holds and perfect knowledge of the bias, the 
measurement significance becomes c 7 /7 ~ 0.015. We 
extrapolate these results to a full-sky galaxy-redshift sur- 
vey, assuming a scale- independent J- = /, to estimate the 
measurement significance of T to about ajr jT ~ 0.05 for 
k < 0.03 h/Mpc, which can be a factor of 6 times stronger 
for fc < 0.1 h/Mpc. 

A similar multi-tracer analysis with the shot-noise can- 
cellation technique but without weak lensing measure- 
ments has been performed [33| to estimate the measure- 
ment significance of the general relativistic coefficients TZ 
and V in galaxy clustering. On large scales, the proper 
relativistic formula in Eq. must be used in galaxy 
clustering and in the power spectrum analysis, Eq. (|92[) 
is conveniently expressed as Eqs. (|93|) and (f94|) with the 
two coefficients in Eqs. (|9"5"|) and In an all-sky galaxy 
survey at z = ~ 1 with fc < 0.03 /i/Mpc, it is found [33[ 
that the galaxy samples with e = 3 and p = 0.4 can con- 
strain the general relativistic effects to 



au/Tl^ 0.1 and a v jV ~ 10.0 



(106) 



if all the halos of M > lO^h" 1 M Q can be utilized. 

Compared to the constraint on ct 7 /7 in the former 
analysis, the constraints on the relativistic effects in the 
latter analysis are weaker, simply because only modes 
close to the horizon scales are sensitive to the relativistic 
effects and also because the former analysis makes use of 
additional information from weak lensing measurements. 
Note that the measurement significance in TZ is about 
half of what is expected for J 7 . However, since the two 
quantities are proportional with an enhancement of TZ 
over T given by Eq. (|105[) that is different for each model 
and roughly of order 3/2 ~ 3, the deviation in TZ may 
be enhanced, and the measurement significance thereof, 
\A7Z/cr-fi\, can become comparable to \AT/ajr\ or even 
exceed it (see Figs. [S][7] and Table IJ). 

To approximately estimate the measurement signif- 
icance of deviations in the relativistic galaxy-density 
fluctuation Eq. (|94|) in non-standard cosmologies from 
its ACDM counterpart, for simplicity, we assume scale- 
independence in AJ 7 , ATZ, and AT 5 , i.e., we adopt the 
sub-horizon limit of the corresponding model at all scales. 
More specifically, we quote the measurement significance 
using predictions at fc ~ 20 Hq and z = 0. We summarize 
our results in Table HI where we fix the cosmological and 
model parameters to the values defined in ^TTl assum- 
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ing that they are accurately determined, for instance, by 
Planck. In this simplified approach, the deviation of up 
to > 20% from the ACDM prediction that we observe in 
the velocity coupling term 1Z in certain quintessence and 
modified gravity models (see Fig. [7]) may be detectable at 
the > 2cr-level, which is comparable to the measurement 
significance expected in the deviations A J 7 (see Tabic [J). 
Thus, the relativistic corrections contain additional in- 
formation on gravity and dark energy, which needs to 
be taken into account in consistent horizon-scale tests of 
departures from ACDM using horizon-scale galaxy clus- 
tering. This conclusion even holds in the case of includ- 
ing more small-scale modes (k < 0.1 h/Mpc) and perfect 
knowledge of bias for constraining J 7 . 

However, we stress again that in our estimation of the 
detection significance, we restrict to a simplified descrip- 
tion of the modified gravity theories and dark energy 
models, neglecting the strong scale dependence of T . 1Z, 
and V at k < 10-ffo, which wc observe in a number of 
modified gravity cosmologies. A proper relativistic de- 
scription within the PPF framework is, therefore, essen- 
tial for constraining those models on horizon scales. We 
leave a more sophisticated analysis of the implications of 
these scale dependencies on the constraints inferred from 
a multi-tracer analysis of future galaxy-redshift surveys 
to future work. In this context, it will be of great in- 
terest to evaluate whether modifications of the nature 
described by, e.g., the PHEN-1 model, which has been 
designed to reproduce ACDM at quasistatic sub-horizon 
scales but shows deviations at near-horizon scales, will 
be distinguishable from the concordance model through 
a horizon-scale test of the kind described here. 



V. CONCLUSION 

We explore the signatures of modified gravity and dark 
energy models in the Newtonian [86| and relativistic con- 
tributions [25l432l | to galaxy clustering. Thereby we adopt 
the Hu-Sawicki PPF formalism Q, which unifies the lin- 
ear perturbations around the FLRW background for a 
large class of modified gravity and dark energy models 
and provides a consistent and efficient framework to de- 
termine the linear perturbations from the near-horizon to 
the quasistatic sub-horizon scales. As example cosmolo- 
gies, we consider a number of quintessence models, f(R) 
gravity, a Brans-Dicke universe with dynamical Brans- 
Dicke parameter and scalar-field potential, both branches 
of the DGP model, and phenomcnological modifications 
of gravity. Since for general scalar-tensor theories, a PPF 
description in the framework of [(| has not been formu- 
lated previously, we develop a calibration procedure for 
finding PPF functions that determine the linear fluctu- 
ations of minimally- and non-minimally-coupled scalar- 
tensor gravity models. We test our approach with a 
number of quintessence and extended quintessence mod- 
els, finding good agreement between the predictions of 
the PPF perturbations and the model-specific fluctua- 



tions. In the limit of f(R) gravity, which is equiva- 
lent to a scalar-tensor model with Brans-Dicke param- 
eter uj = 0, we find that our fitting prescription repro- 
duces the known f(R) PPF fits of @. Our PPF cali- 
bration fits for general scalar-tensor theories provide an 
efficient procedure for consistently computing the asso- 
ciated perturbations on all linear scales. They improve 
upon the widely-used simple quasistatic sub-horizon ap- 
proximation that can be applied to obtain perturbations 
of scalar-tensor models at large scales, which, however, 
breaks down near the horizon and cannot correctly de- 
scribe the clustering of the dark energy associated with 
the scalar field of the minimally-coupled models on large 
scales and its additional horizon-scale anisotropic stress 
in the non-minimally-coupled cases. 

We have applied the PPF formalism to galaxy clus- 
tering, providing a proper relativistic treatment on hori- 
zon scales, where the relativistic effects become substan- 
tial. As in GR, the relativistic formula of galaxy clus- 
tering in modified gravity theories can be characterized 
by three parameters: the Newtonian growth rate J 7 , 
the velocity contribution 7Z, and the potential contribu- 
tion V . In modified gravity, the Newtonian growth rate 
T = — knVm/ A m slightly differs from the usual logarith- 
mic growth rate of structure / and the redshift-space dis- 
tortion in galaxy clustering measures J 7 , not the usual /. 
Compared to GR, in modified gravity and dark energy 
models, the three quantities J-, 7Z, and V have different 
values and exhibit scale dependencies. 

We approximately estimate the detection significance 
of these effects in a multi-tracer analysis of future all- 
sky galaxy-redshift survey data employing halos of mass 
M > 10 10 /i _1 Mq at low rcdshifts. Such an analysis 
may have the potential of discriminating between the 
relativistic correction sourced by the velocity coupling 
produced in a number of quintessence or modified grav- 
ity scenarios and their counterpart from the concordance 
model. Thereby, we find a significance level which is 
comparable to the measurement significance expected in 
the deviations of the logarithmic matter-density growth 
rate. Thus, we conclude that the relativistic corrections 
contain additional information on gravity and dark en- 
ergy, which needs to be taken into account in consistent 
horizon-scale tests of departures from the concordance 
model using horizon-scale galaxy clustering. 

Note that we restrict to a simplified estimation of the 
measurement significance, neglecting strong scale depen- 
dencies in the velocity fluctuation at the largest mea- 
surable scales, which appear in a number of modified 
gravity cosmologies. In this context, we further design 
modifications of the concordance model, which repro- 
duce the standard relations at quasistatic sub-horizon 
scales and cannot be tested by conventional cosmological 
probes. The multi-tracer analysis, however, may pro- 
vide a method to discriminate between those scenarios. 
We leave a more sophisticated analysis of the parameter 
constraints that may be obtained from future galaxy sur- 
veys to future work. Such an analysis will include effects 
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from additional scale dependencies in the growth and 
in the relativistic coefficients as well as an evaluation of 
the employment of these effects to discriminate between 
horizon-scale deviations and the concordance model. 
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Appendix A: Numerical computation of 
perturbations 



In the following, we shall give details on the numer- 
ical integration of the PPF perturbations. In order to 
calibrate the linear PPF perturbations for scalar-tensor 
theories and for the comparison of the performance of 
their corresponding fits, we also need to determine the 
full model-specific fluctuations. For this purpose, we 
briefly review the linear perturbation theory of scalar- 
tensor gravity models of the kind defined by Eq. $Z§. We 
refer to (37j for a more comprehensive discussion of this 
modified cosmological perturbation theory. For simplic- 
ity, we restrict the linear perturbations to a spatially-flat 
matter-only universe. 



2. Minimally-coupled scalar-tensor gravity 
perturbations 

After some algebraic manipulation, the perturbed 
Einstein equations, scalar-field equation, and energy- 
momentum conservation equations of the minimally- 
coupled scalar-tensor (or quintessence) models in a 
spatially-flat FLRW background (see mi A 1[) become 

-$ - ^ttkPt -z-St, 

k H 2 



vL 

Sip' 



2H 2 
2/T 2 

-k H § - v m , 



-k 2 H 



H 



2fc| 
V>' 



(A4) 
(A5) 

Given the background from £]II A 11 i.e., H or w e s, U, 
and <p' , the evolution of these perturbation variables are 
solved by numerical integration over In a with the initial 
conditions at early time a\ < 0.01 corresponding to the 
matter-dominated era in GR, 

= (A6) 



-p'Sp, 



(A2) 
(A3) 



A, 



2 



Spi = o, 
* - f 6, 



(A7) 
(A8) 
(A9) 



where the comoving-gaugc curvature Q is a constant. 



3. Non-minimally-coupled scalar-tensor gravity 
perturbations 



1. PPF perturbations 

Given the background expansion history H(z) or 
w e g(z) and the PPF functions and parameter, g, fa, fe, 
and cr, we evolve the linear perturbations in the PPF for- 
malism described in through the conservation equa- 
tions, Eqs. ((56)) and (|57)) . the differential equation for V, 
Eq. (|63p , with source Eq. (|60p , and the modified Poisson 
equation, Eq. (|52"j) . They compose a system of coupled 
first-order differential equations in A m , V m , and T as well 
as an additional constraint equation for The initial 
conditions are set at early times ai < 0.01 with 

V mti = \k H ^j [(g> + + (1 - g + .9')$-] ) , 

(Al) 

where the initial metric perturbation is a constant 
and , = T = 0. 



For the non-minimally-coupled scalar-tensor (or gen- 
eralized/extended quintessence) models in Brans-Dickc 
representation of Eq. ([2]) as well as for a perturbed 
spatially-flat FLRW background metric, the perturbed 
parts of the modified Einstein equations arc (cf. (37j ) 

Sip 

H 2 pV m 
H 2 k H 



$ + * 



(A10) 
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From the perturbed scalar-field equation, one obtains the 
dilaton fluctuation equation 




Finally, the conservation equations imply 



k H 



y m 

k H 



-k H V m , 
V'+V m = k H t>. 



(A14) 
(A15) 



To solve for the non-minimally-coupled scalar-tensor 
model perturbations, given the background from mi A 21 
i.e., H or w c {[, tp, and U, we combine the perturbed mod- 
ified Einstein and energy-momentum conservation equa- 
tions to obtain a second-order differential equation for 

which includes dependencies on Sip. We combine this 
with the dilaton equation depending on <£> and Sip. We 
use the perturbations for general F, U, and Z, with- 
out rewriting it in the Brans-Dicke representation (see, 
e.g., We follow [53| in this combination, how- 

ever, solving the system of coupled differential equations 
with respect to In a instead of cosmic time. This pro- 
duces two coupled second-order differential equations for 
$ and Sip, which we solve with the initial conditions at 
early times a[ < 0.01 corresponding to the GR limit, i.e., 
$| = 5tp[ = 5ip[ = 0, where the initial metric perturba- 
tion $; is a constant. 

Note that since we use divisions by Z, we cannot set 
Z = to strictly recover f(R) gravity. For f(R) gravity, 
we therefore solve the perturbations given in [54| for the 
designer model with a code developed and tested in (5|| . 
Furthermore, we use a code, where we solve the non- 



minimally-couplcd scalar-tensor perturbations given here 
in Brans-Dicke representation setting u = 0. We test the 
three non-minimally-coupled scalar-tensor perturbation 
theory solvers in the limit were uj — > and find good 
agreement between the codes. 



Appendix B: Numerical computation of the NCBD 
background 

Finally, note that in order to obtain the background 
for the NCBD model, we follow (53[, assuming matter 
dominance and using the first Fricdmann equation, the 
scalar- field equation with dH/dt replaced by the second 
Fricdmann equation, and the energy conservation equa- 
tion, i.e., Eqs. (fT9|) through (|22|) . to solve for <p, a, and 
p m as functions of cosmic time t. This produces coupled 
differential equations that are of second order in ip and 
first order in a as well as p m , which we integrate over t 
using the initial conditions 



t 



K Pm,i 



4 



(Bl) 
(B2) 



with constant initial scale factor a; and constant ipi as 
well as pi for the scalar field, where the dot represents a 
cosmic-time derivative. 
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